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Abstract 

We consider a fixed number of noninteracting bosons in a harmonic trap. 
The determination of the exact microcanonical ground state number fluctua- 
tion is a difficult enterprise. There have been several theoretical attempts to 
solve this problem approximately, especially in ID models where analytic re- 
sults were found using some asymptotic formulae from number theory. Here, 
we obtain the exact number fluctuation curves, and show that these exact 
curves are substantially different from the results presented in literature for 
small particle number and finite temperature. Further, we show how these 
approximate microcanonical results may be derived using a canonical ther- 
modynamics approach. 
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I. INTRODUCTION 



As is well known, the traditional approach of determining the fluctuation of single- 
particle occupation numbers fails for bosons as the temperature T — > 0. It predicts that the 
fluctuation tends to infinity in this limit, rather than zero as described by a real system. 
After the experimental success in producing Bose-Einstein condensation of trapped alkali- 
metal atoms, there is renewed interest in calculating the number fluctuation using either 
canonical or microcanonical ensemble approach Most of these, however, use some 

forms of approximation to the real microcanonical fluctuation. In 0] the authors define the 
so-called Maxwell's Demon ensemble, and calculate the microcanonical fluctuation for a 3D 
trap numerically. In |7| we formulated a method of calculating the exact microcanonical 
fluctuation using combinatorics for a perfectly isolated BEC with finite number of particles 
which, to our best knowledge, has not been presented before. 

In this report we focus on a system of bosons in an one-dimensional harmonic trap. The 
number fluctuation of this system was calculated by the authors in |L| , using the well-known 
asymptotic formulae from number theory to approximate the exact expansion coefficients 
of the N-particle partition function. Their method yields results which were supposed to 
be microcanonical. However, we shall show in this brief report that their approximate 
method yields the same results one obtains from canonical ensemble averaging, or a related 
formulation by Parvan et al [11|. At low temperatures and for finite number of particles (N < 
100), these results are very different from the ones calculated using the exact combinatorics 
method. 

The layout of this paper is as follows. In section H we establish the notations and 
outline the method of calculating the exact ground state number fluctuation Section 
|H contains the main focus of this brief report in which we derive the results of |Q using 
thermodynamics, and show that they agree with the ones obtained from canonical method 
as given by ||11|| . Finally, we compare the exact combinatorics results from section [II] and 
the canonical results from section [m[ 



II. MICROCANONICAL FORMALISM 

We consider an isolated system with discreet energy levels consisting of N particles in the 
ground state at zero temperature. Given an excitation energy E, there can be one particle 
which absorbs all the energy and gets excited to an excited state, or two particles which share 
the available energy... etc. The number of excited particles N ex can be somewhere between 
1 and N, depending on the excitation quanta, such that the total number of particles N 
remain constant. This defines the microcanonical ground state number fluctuation. Denote 
co(E, N ex , N) to be the number of possible ways of distributing E among exactly N ex particles, 
then the probability of exciting exactly N ex particles in a N-particle system at excitation 
energy E is given by f7|,§]: 

P(E,N a ,N)= (!) 

Y. N N „=ME,N„,N) y 

Expression (|I|) is normalized and has the following properties: 
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P(0,N ex ,N) = 5 0Nex , (2) 
P(E,N ex ,N) = N ex > N. (3) 

The moments of the probability distribution now read: 



N 



(N ex ) = ]T N ex P(E,N ex ,N) (4) 

N ex =l 
N 

N 2 ex ) = E N 2 ex P(E,N ex ,N) (5) 



N ex =l 

The number of particles in the ground state Nq and the excited states N ex are related via 
N = Nq + N ex . Nq and N ex are allowed to vary while N is fixed so that the number 
fluctuation from the ground state is given by: 

((5N f) = (N 2 ex ) - (N ex f 

> 2 ) " (A^o) 2 (6) 



Expression (|]) is true in general regardless of the ensemble and statistics of the particles. 

Next, we specialize to a one-dimensional harmonic trap. Let us put hu = 1, then the 
excitation energy from the ground state is E — n. The partition function of the system is 
well-known 

Zi(P) = ^— (7) 
1 — x 

Zn(P) = U T ^- (8) 

where x = . This is just the generating function pertaining to the restricted linear 
partitions of an integer n in number theory ||. Expressing (§) in power series of x 

CO 

Z N (f3) = Y,n(n,N)x n (9) 

n=0 

The coefficients Q(n, N) is the number of partitions of n which have at most N summands. 
For example, if n = 5 and N = 4, then Q(5, 4) = 6 since 5 = 5, 1 + 4, 2 + 3, 1 + 1 + 3,1 + 
2 + 2,1 + 1 + 1 + 2. In connection to our system, Q(n, N) is interpreted as the number of 
possible ways to distribute n quanta among utmost iV particles. Clearly, 

N 

tt{n,N) = £ u{E,N ex ,N) (10) 

N ex =l 

While the coefficients Q(n, N) may be found analytically by expanding Z^(/3) in series, 
the exact determination of the microstate lo(E, N ex , N) is a difficult problem. We may 
find uj(E, N ex , N) using combinatorics as described in detail in 0, or, alternatively, use the 
well-known identity || 

u(n, N eX) N) = n(n, N ex ) - n(n, N ex - 1) (11) 

Using (|TT|), (|l]) and the exact microcanonical ground state fluctuation may be deter- 
mined. 
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III. CANONICAL FORMALISM 



The ground state number fluctuation may in general be expressed in terms of the mo- 
ments of occupation numbers as [ID),!!]!!: 

W) = E " (n*> 2 ) (12) 

k 

The sum runs through all the allowed k values in the ground state defined at zero tempera- 
ture. For bosons, only k = state applies. Canonically, the moments are given by: 

1 N 



("*> =— Xy^-iCfl (13) 

1 



Z N 7= i 



r- 

N 

ni) —Y, [2j'-l] Z v ,( i) (14) 



Z N 



We can also obtain the canonical fluctuation using a thermodynamics approach . Given 
the N-particle partition function (|8|), the entropy of the system reads 

N ■ N 

i=l I- 1 J i=l 

iV 

= (3E - E ln[l - x j ] (15) 

Note that the canonical multiplicity Q(/3, N) by definition is given by: 

Q((3,N)=exp[S N ((3)} =x~ E Z N {f3) (16) 

Compare ( |16"D with (|9|), we see that the canonical multiplicity is a one-term approximation to 
the exact value of the expansion coefficients of the series @. To make a direct comparison 
with the results given in Q, we approximate the sums in (|l^). The first term gives the 
energy-temperature relationship: 

E = H = § [x-/- 1] ^ /o [x-/- 1] 

To retain the N-dependence we use the Euler-Maclaurin sum formula in approximating the 
second term: 

N 7T 2 1 1 

Em-*^- y - + -^ as) 

Using (|ll|) and fll5D , the canonical probability of finding N ex particles in the excited states 
is given by: 
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P(P,N ex ,N) 



n(P,N ex )-n((3,N ex -i) 

1 dQ(f3,N ex ) 

n(p t n) dN ex 

1 d exp[S(fJ,N ex )} 
exp[S((3,N)} dN ex 



e P Nex exp 



I ( e -f>N„ _ e -m 







(19) 



Using fll7| ) and the transformations 



;7T, Xr 



c r 



2v^ 



ln(y/ri) 



then, 



P(n,N ex ,N) 



I exp 


- 2 -e 

c 


XN ex _ 




exp 


-l e -XN 

c 





(20) 



This is identical to expression (14) given in ref. [JTJ] which the authors claim give the mi- 
crocanonical fluctuation. That the canonical method given by (|12|), (|l3l) , and (0) and 
the above method are equivalent is shown in Fig. (|1|) J7|. Both yield the same canonical 
fluctuations. 



IV. MICROCANONICAL VS. CANONICAL 

Fig. (0) compares the relative fluctuations for N = 10, 50, 100 particles. The respective 
curves of both ensembles tend to agree as iV and n get large (ie. thermodynamic limit), 
otherwise the results are very different. 

It is clear from the discussion above that the method employed by the authors in ref. 
yields canonical results. Comparison between their method and the thermodynamics method 
makes it clear on why this is so: both use what essentially is a saddle point approximation to 
the exact multiplicity Q(n, N) (p7v(m) in their notation) from the series <^). Only when the 
multiplicity found exactly are the results microcanonical. This clearly makes it more difficult 
to incorporate interaction, and treat the true condensate fluctuations of weakly interacting 
isolated bosons microcanonically. 

The author thank R. K. Bhaduri and M. V. N. Murthy for fruitful discussions. This work 
was supported by NSERC grant # 199160-01 and NSERC scholarship # PGSA-221708-1999. 
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FIGURES 



FIG. 1. Canonical ground state number fluctuation for N = 1000 bosons in a ID harmonic 
confinement. The solid line shows the fluctuation using formulae ©, (|l|), and (|14|). The diamonds 
were obtained using the thermodynamics method as described in section |III] . 



FIG. 2. Comparison between canonical and microcanonical relative fluctuations. The critical 
temperature is given by T c = N/lnN with fiuj = 1 [1]. The inset shows a close-up of the low 
temperature part. 
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